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ABSTRACT 

O 

; There has long been evidence that low-mass galaxies are systematically larger 

in radius, of lower central stellar mass density, and of lower central phase-space 

a ■ 

density, than are star clusters of the same luminosity The larger radius, at a 
comparable value of central velocity dispersion, implies a larger mass at similar 
luminosity, and hence significant dark matter, in dwarf galaxies, compared to no 
dark matter in star clusters. We present a synthesis of recent photometric and 
>■ ■ kinematic data for several of the most dark-matter dominated galaxies. There is 

a bimodal distribution in half-light radii, with stable star clusters always being 
smaller than ~ 30pc, while stable galaxies are always larger than ~ 120pc. We 
extend the previously known observational relationships and interpret them in 
terms of a more fundamental pair of intrinsic properties of dark matter itself: 
dark matter forms cored mass distributions, with a core scale length of greater 
than about lOOpc, and always has a maximum central mass density with a narrow 
range. The dark matter in dSph galaxies appears to be clustered such that there 
is a mean volume mass density within the stellar distribution which has the very 
low value of about 0.1 M & pc~ 3 (about 5GeV/c 2 cm -3 ). All dSphs have velocity 
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dispersions equivalent to circular velocities at the edge of their light distributions 
of ~ 15km s _1 . In two dSphs there is evidence that the density profile is shallow 
(cored) in the inner regions, and so far none of the dSphs display kinematics 
which require the presence of an inner cusp. The maximum central dark matter 
density derived is model dependent, but is likely to have a mean value (averaged 
over a volume of radius lOpc) of ~ 0.1 M pc~ 3 (about 5GeV/c 2 cm -3 ) for our 
proposed cored dark mass distributions (where it is similar to the mean value), 
or ~ 60 M Q pc~ 3 (about 2TeV/c 2 cm -3 ) if the dark matter density distribution 
is cusped. Galaxies are embedded in dark matter halos with these properties; 
smaller systems containing dark matter are not observed. These values provide 
new information into the nature of the dominant form of dark matter. 

Subject headings: dark matter — galaxies: individual (dSph) — galaxies: kine- 
matics and dynamics — Local Group — stellar dynamics 



1. Introduction 

The distributions of total luminosity and of central stellar velocity dispersion for star 
clusters and for dwarf galaxies overlap, so that the faintest galaxies have approximately the 
same values of these physical parameters as do star clusters, with galaxy luminosities ex- 
tending as faint as ~ 1O 3 L , with line-of-sight central velocity dispersions of ~ 10 km/s. The 
half-light radii (radius containing one-half the total luminosity) of the galaxies, however, are 
significantly larger (hundreds of parsec) than those of star clusters (at most tens of parsec). 
This leads, through the virial theorem, to significantly larger inferred masses for the dwarf 
galaxies, compared to star clusters of the same luminosity and velocity dispersion. Indeed, 
the derived values of central and global mass-to-light ratios for the gas-poor, low-luminosity, 
low-surface brightness satellite galaxies (classified as Dwarf Spheroidal galaxies; dSph) of the 
Milky Way are high, up to several hundred in solar units, making these systems the most 
dark-matter-dominated galaxies in the local Universe (e.g. Mateo 1998 for a convenient re- 
view of early work). As we discuss further below, they are th e ideal test-beds for constrain ing 



the nature of the dark matter that dominates their gravity (jOstriker fc Steinhardtll2003l ). 



The dSph galaxies and star clusters share a further observed property, that organised 
orbital rotational energy of the member stars is negligible compared to the energy in dis- 
ordered motion, which is measured by the stellar velocity dispersion at a given location. 
Similarly to pressure gradients in a fluid, the stellar velocity dispersion provides the support 
against self-gravity, but unlike the fluid case, stellar pressure can be anisotropic, generating 
galaxy shapes which need not be spherical. Systems in which angular momentum support 
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against gravitational potential gradients can be ignored when analysing the kinematics of 
member stars are designated as 'hot'. 

It has been known for the past twenty years that there are well-defined, and probably 
fundamental, scaling relations between the half-light radius (or core radius), the central 
velocity dispersion and the luminosity of hot stellar systems (e.g. Kormendy 1985; Bender, 
Burstein & Faber 1992; Zaritsky, Gonzalez & Zabludoff 2006a,b). It has further been long- 
established that the globular star clusters in the halo of our Galaxy show distinctly different 
scalings from the dSph galaxies, and that the dSph galaxies in turn have different scalings 
from more luminous hot galaxies (e.g. Kormendy 1985, his Figure 3; Burstein et al. 1997). 
Dynamical effects over their long lives have modified the size and luminosity distributions 
of the Galactic globular clusters (e.g. Fall & Rees 1977; Gnedin & Ostriker 1997), so it 
is important that robust studies include star clusters of all ages and in all environments, 
including globular star clusters in external galaxies, nuclear star clusters, and young massive 
star clusters, significantly younger than globular clusters (e.g. Walcher et al. 2005; Seth et 
al. 2006). 

The specific combination of central velocity dispersion (o"o) and half light radius (r^) 
r^dQ 1 oc PhC 3 , is a convenient measure of the phase-space density, where ph is the mean 
density within a half-light radius. The similarity of velocity dispersion for star clusters and 
dwarf galaxies of the same luminosity, combined with the factor of ~ ten difference in their 
half-light radii, implies a systematic difference of some two orders of magnitude in the value 
of the phase-space density at fixed stellar mass (see e.g. Walcher et al. 2006). 

What is the physical explanation for these differences between star clusters and low- 
luminosity galaxies? Clearly, the presence of dark matter in dSph galaxies, but not in star 
clusters, is a critical distinction, and provides an opportunity to identify underlying physics 
of Dark Matter. For example, the suggestion by Mateo et al. (1993) that there is an apparent 
minimum dark halo mass of ~ 2 x 10 7 M Q , deduced from the available dynamical studies, 
implies a small-scale limit to the dark matter power spectrum unlike that assumed in ACDM 
models (e.g. Moore et al. 1999; Klypin et al. 1999). The minimum halo mass suggestion was 
shown to still be valid in a significantly extended sample, including dSph satellites of M31, in 
addition to those of the Milky Way, in an important study by Cote et al. (1999; see especially 
their figure 3). Cote et al. also provided one of the earliest robust demonstrations that the 
internal kinematics of dSph galaxies are in general unaffected by external tidal forces from 
their host galaxies, so that the results from application of equilibrium dynamical analyses 
are reliable. The putative minimum dark-halo mass was still found to be a pprop riate, in 
the larger sample wit h more extensive data reviewed by IWilkinson et al. I (120061 ) and by 



Gilmore et al. I (120061 ) 
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Within the broad class of dwarf galaxies, of which dSph are the least luminous mem- 
bers, one can apply simple models to the data and obtain scaling relations between such 
quantities as derived central dark matter density and observed central velocity dispersion 
(e.g. Kormendy & Freeman 2004). The dSph galaxies are systematically discrepant in their 
correlation fit, falling below the extrapolated trend to larger central densities as luminosities 
decrease. We provide an explanation here by showing that the dSph form the limit of such 
relations, not a continuation. These correlations can be used to consider compatability with 
various parameterizations of the power spectrum of primordial density fluctuations. Dwarf 
galaxies play a special role, in that they appear to be the smallest systems in which dark 
matter dominates, and so provide a powerful test of the power spectrum on the smallest 
scales. The smallest scales on which dark matter particles cluster depends on the physical 
characteristics of the dark matter itself (e.g. Green, Hofman & Schwarz 2005). Determining 
this smallest scale is the goal of the present analysis. 

In this paper we revisit the established correlations and scaling relations for dwarf 
galaxies and for star clusters. Stellar velocity data now exist for stars across the face of 
several of the dSph, allowing an analysis that goes well beyond that possible with just the 
central value of the velocity dispersion, a limitation in early studies. The discussion below 
takes account of these new data, where available. Recent imaging data allow a re-evaluation 
of the sizes of star clusters and galaxies, strengthening the case for a real discontinuity 
between star clusters and galaxies. We interpret our findings in terms of a more fundamental 
pair of intrinsic properties of dark matter itself. 



2. The sizes and internal kinematics of star clusters and galaxies 

The existence of a clear observational distinction between massive star clusters and low- 
mass galaxies has been substantially strengthened recently, both through detailed studies 
of more luminous and massive star clusters in a wide range of environments, and through 
discovery of a large number of extremely low-luminosity satellite galaxies around the Milky 
Way (e.g. Willman et al. 2005b; Belokurov et al. 2007; Zucker et al 2006a, 2006c) and around 
M31 (e.g. Zucker et al 2004, 2006b), mostly based on imaging from the Sloan Digital Sky 
Survey (SDSS). Fig. 1 shows the current sample in a plot of half-light radius against absolute 
magnitude in the V-band. Star clusters in all studied environments, with luminosities over 
the whole range from My = —4, L~ 10 3 L , up to My = —15, L~ 1O 9 L , and a wide 
range of ages, invariably have characteristic scale sizes r\ less than about 30pc. Thus the 
dynamical range over which both star clusters and galaxies exist, and over which there is a 
distinct size dichotomy, has now been established to cover some six orders of magnitude in 
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stellar luminosity. Available direct studies of the stellar initial mass function in star clusters 
and dSph galaxies (e.g. Wyse et al. 2002) show this range in luminosity corresponds to a 
similar dynamic range in baryonic mass. 

Figure 1 makes evident that there is a robust maximum radius to star clusters, at all 
luminosities. Figure 1 also illustrates that the dSph galaxies in both the Milky Way and 
in M31 have a minimum characteristic radius, and this minimum is a factor of > 4 larger 
than the largest star clusters. With the exception of the recently discovered object ComaBer 
(Belokurov et al. 2007) - the largest and brightest of the three systems indicated by ringed 
red circles in Fig. 1 - which manifestly merits further study, but is one of the two new dSph 
which lie between the Sgr dwarf and the Magellanic Clouds, and which show significant 
indications of tidal disruption, there is no known object in the size-gap between ~ 30pc and 
~ 120pc. 

It is more correct to say that, modulo ComBer, there is no known stable object in the 
size gap. This intermediate size might be occupied transiently by a larger object (a dwarf 
galaxy) in the very late stages of disruption by external (Galactic) tides, or by a small object 
(globular cluster) in the last stages of evaporation. In the first of these low velocity 

dispersion compact core can be generated transiently, if outer, hotter, stars are removed by 
a suitable tide, while in the second case the density profile changes systematically from the 
small value typical of a compact star cluster to a a very large value, almost constant density, 
covering all possible radii during that (short-lived) process. 



Willman 1 (radius 20pc; IWillman et al. 1120061 ) and Segue 1 (radius 30pc; lBelokurov et al. 



20071 ) are also newly discovered and interesting tests of the conclusions of this section, with 
Segue 1 showing evidence for significant tidal disruption. The two largest Galactic globu- 
lar clusters are Pal 14, with size 28 pc, and Pal 5, 24 pc, w hich is in an advanc e d sta ge 
of tidal disruption with prominent streams of stripped stars ( lOdenkirchen et al. I (120031 )). 
The largest Ultra Compact Dwarf galaxies (UCDs), also with size 25pc, are associated with 
the centre of the Virgo cluster, and the galax y M87, and were for a long time suspected 
(IDrinkwater et al.l 120031 ; lHasegan et al. 1 120051 ) of being small galaxies severely affected by 
tides. Our interpretation of their status, based on their position in Figure 1, is that they are 
simply very massive star clusters, with no associated dark matter. 

There are two ve ry recent detailed dynami cal analyses of the masses of ultra com- 
pact d w arf galaxies, bvlHilker etal I (120071 ) and by lEvstigneeva. Gregg. Drinkwater. fc Hilker 
(120071 ). iHilker etal I (120071 ) de rived dynamical masses for five ultra-comp a ct dw arfs and one 
dE nucleus in Fornax, while lEvstigneeva. Gregg. Drinkwater. fc Hilker! (120071 ) studied six 
Virgo UCDs and five very luminous Fornax UCDs. They show that all these systems are 
similar, in structure and dynamics, and that the dynamical mass-to-light ratios for the UCDs 
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are consistent with simple stellar models: there is no evidence for any dark matter associated 
with these stellar clusters. They are the (very) high- mass/high-luminosity extreme of more 
typical globular cluster populations. 

There are systematic effects that need to be taken into consideration when interpreting 
Figure 1. First, the half-light radius (whose definition is the obvious one, the radius enclosing 
one-half the total luminosity) can be robustly estimated only in systems with a well-defined 
and convergent luminosity profile. In many cases the parameter published is a 'core radius', 
that radius at which the projected surface brightness has fallen to one-half its central value. 
For the commonly used King-model fits to star clusters, the core and half-light radii are 
similar, except in (low concentration) cases where the object has an extended tail to the 
brightness distribution. In this case, the derived core radius underestimates the true half- 
light radius. For a modelled projected surface brightness fi(r) described by 

Mr)=/, (l + r 2 /a 2 )" 7/2 , (2-1) 

where a is the scalelength of the core, /xo the central surface density and 7 the power-law 
decline of the surface density at large radii, King models have 7 ~ 2, and the Plummer 
sphere has 7 = 4. The projected integrated luminosity L p is 



27r/i 

where the integral to infinity is 



MO = ^(a 2 -o>(° 2 + <T^). (2-2) 



L M = ^L. (2-3) 

7-2 

For the Plummer sphere L p (a) = l/2L tot . That is, the physical meaning of the scale param- 
eter a in this Plummer case is a half-light radius. More generally, the relation between the 
scale parameter a and the core radius r c , defined as above, can be shown to be 

r c = a(2 2 ^ - lf'\ (2-4) 

where a is the (cylindrical) radius which encloses one-half the total luminosity, so that for a 
Plummer sphere observed in projection the half-light radius a and core radius r c are related 
by a ~ 3r c /2. Core radii, where fitted, are adequate approximations to half-light radii for 
King model globular star clusters, and are lower limits to half-light radii for dSph galaxies, 
while the Plummer scale parameter is a half-light radius, to the accuracy of a Plummer 
model fit to the data. Conservatively, we adopt and identify three cases of core radii for 
dSph galaxies at face value, rather than converting to larger, but more uncertain, half-light 
radii. 
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Second, only the extent of the baryonic component is being measured, and there is no 
guarantee - or need - for mass to follow light. Given that most, if not all (as we argue 
here), galaxies are embedded in extended dark matter halos, photometric determinations 
are probably a lower limit on the scale length of the total mass distribution (we discuss this 
below for two specific caes, UMi and Fornax). Figure 1 is therefore probably conservative in 
describing mass: an even larger distinction between star clusters and dSph galaxies would 
be seen were one able to plot parameters describing mass rather than light. 

Figure 1 is further conservative, in that it is possible that the recently discovered very- 
low luminosity dSph are in fact larger than is shown. The recent history of observational 
studies of nearby low- luminosity galaxies, in which individual stars are resolved and the 
extent on the sky is measured through star counts, corrected for foreground stars in the 
Milky Way Galaxy, has been that their radial extents (a nd hence total luminositie s and 



half-light radii) tend to be somewhat under-estimated (eg lOdenkirchen et al. I (120011 ) ) . As 
photometric data are extended to lower surface brightnesses, and as kinematic studies of 
individual member stars develop, allowing foreground stars to be rejected on the basis of 
line-of-sight velocity, the galaxies are typically found to be larger than first measured. For 
example, for the best-studied dSph the estimated total extent has changed by a factor 1.5-2 
over the last 10 years. In contrast, star clusters really do have steep outer light profiles, 
and thus have well-defined observational parameters. Both these effects lead to a systematic 
observational underestimate of any gap in spatial scale between star clusters and galaxies. 



3. Masses and mass distributions, cores and cusps 



It has long been known that the observed value of ~ 10 km/ s for the line-of-sight velocity 
dispersions of Local Group dwarf spheroidals (dSphs), together with their ~ 200pc half-light 
radii, implies mass-to-light ratios M/L of up to ^ 1OOM /L . Until recently, most of these 
estimates of Mj L were based upon a measurement of only the central value of the velocity dis- 
persion, and upon the assumption that the mass profile follows the light profile. The availabil- 
ity of datasets of radial velocities for hundreds of individual stars spread out in radius across 
the nearby dSphs, obtained by several groups, has changed all this. To date, the line-of-sight 
velocity dispersion profiles in the Fornax, Draco, Ursa Minor, Carina, Leo I, Leo II, Sextans 
and Scul ptor dSphs have been mapped to the (rather poorly defined concept of an) opti- 



cal edge (IMateo 



Muhoz et al. 



19971: iKlevna et al II2OO2I: 



2005: 



Sohn et al. 



20061 : 



Wilkinson et al. 



Koch et al. 



2004 



Wi 



20061 : iBattaglia et al 



kinson et al. 



200 



6; 



2006: 



Walker et al. 



20061 ) - see Figure 2 for a sample of profiles obtained by our group. Central velocity dis- 
persions have been derived for several of the newly discovered extremely low-luminosity 
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dSph satellites of both the Milky Way and M31 [Kleyna et al. 



(2005 


); 


Ibata et al. 


(2006) 



Muhoz et al. I (l2006bl )]. Several major studies are underway so that all the known dSph 
galaxy satellites with —13 < My ^ —8 will soon have much improved determinations of 
their dynamical mass distributions inside their optical radii, while some information will be 
available on the several dSph galaxies with —8 ^ My ^ —4. Our conclusions in this paper 
lead us to predict extended dark matter halos for the systems with characteristic radii that 
place them in the 'galaxy' regime of Fig. 1. Future more extensive studies of these very faint 
galaxies, as well as continuing discoveries, will allow our predictions to be tested. 



3.1. Mass modelling for hot stellar systems 

Stellar systems with no net internal angular momentum maintain their scale size by 
the pressure support of random stellar motions against the gravitational potential gradient. 
This pressure is naturally tri-axial in a collisionless system. A robust dynamical analysis 
of such a system involves solution of the Collisionless Boltzmann Equation (CBE) for some 
appropriate (stellar) tracer particle phase-space distribution function, determining the mass 
distribution which generates the gravitational potenti al gradients along wh ic h the stars orbit. 



Such models are being applie d and further developed jKleyna et al. I (120011 ); IWilkinson et al. 



( 120021 ); I Walker et al. I (120061 )] to the few best-studied dSph galaxies. Such analyses are 
appropriate for, and require, velocity and position data for (at least) several hundred tracer 
stars distributed across the dSph. For most dSph galaxies studied to date, the more limited 
data available justifies analysis only of the first velocity moment of the distribution function, 
the velocity dispersion as a function of radius. 

In a collisionless equilibrium system the Jeans' equations are the relation between the 
kinematics of the tracer stellar population and the underlying (stellar plus dark) mass dis- 
tribution. In terms of the intrinsic quantities, and assuming spherical symmetry, the mass 
profile may be derived as: 

r 2 ( 1 Aval n P°t 



M(r) = - — — ^ + 2 



G \v dr r 

where a r (r) is the one- dimensional stellar velocity-dispersion component radially towards 
the center of the mass distribution, (3{r) quantifies the stress term associated with (possibly 
radially variable) velocity orbital anisotropy, and u(r) is the stellar density distribution. 

The quantities directly observed are the line-of-sight velocity dispersion as a function of 
projected radius, R, a(R) =< v(R) 2 > x / 2 and the surface brightness profile as a function of 
projected radius. Given finite amounts of data defining the projected surface brightness and 
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kinematic distribution functions, we can proceed in either of two ways. We may assume a 
priori a parameterised mass model M(r) and velocity anisotropy (3(r) and fit the observed 
velocity dispersion profile; or we may use the Jeans' equations to determine the mass profile 
from the projected velocity dispersion profile, utilising some (differentiable) functional fit to 
the observed light distribution and a (range of) assumed form(s) for the anisotropy f3(r). 
Assuming spherical symmetry, it is straightforward to obtain (of) from the observed line-of- 
sight velocity dispersion using Abel integrals. In what follows, we take the second approach 
to the Jeans' equation analysis: both the spatially binned dispersion profile and the surface 
brightness distribution are fit by an appropriate smooth function, and we assume an isotropic 
velocity dispersion. Figure [3] shows some examples of the fits to the light and dispersion 
profiles used in the analysis. 

It is obvious from the Jeans' equation that radially variable velocity dispersion anisotropy 
is degenerate with mass, making any deductions as to whether or not the inner mass profile 
is cored or cusped in general mo del- dependent. Further information is needed to break this 
degeneracy, and fortunately is sometimes available, as we discuss below. In general how- 
ever, full multi-component distribution-function models using adequately large datasets, as 
discussed in Section 3.4 below, are required to use the information in the data to break this 
degeneracy. 

3.2. Moment equation analyses of inner dark mass distributions 

Jeans' equation dynamical analyses generate three quantities. The most robust is the 
mean dark matter mass density inside the radius where adequate kinematic data are avail- 
able. Similarly robust is the total mass, again inside the radius where adequate kinematic 
data are available. The analysis can also constrain the mass density in a small central re- 
gion, though the limited central spatial sampling makes this less robust, and dependent on 
an adopted underlying mass density profile. It remains the case, of course, that a suitably 
small central cusp in the mass distribution would not be detected so long as it was unresolved 
inside the available kinematic data: very small cusps might be present. 

Our Jeans' analysis mass models are presented in Figure 4. Some caution is required in 
the interpretation of these profiles, given the simplifying assumptions which have been made 
(spherical symmetry; velocity isotropy; smooth dispersion and light profiles). First, given 
that satisfying the Jeans equations is a necessary, but not sufficient, condition for a solution 
of the CBE to be everywhere non- negative (and hence a viable distribution function), the 
models presented here are not guaranteed to correspond to physical models. However, we 
note that for a tracer distribution with an isotropic velocity distribution and density profile 
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v ~ r 7 , the logarithmic slope of any external power- law pote ntial -0 ~ r 15 must satisfy 



5 < 2"f to ensure the non- negativity of the distribution function (lAn Sz Evand 120061 ) . Thus, 
it is reasonable to assume that for a tracer distribution with 7 ~ 0, a cored mass density 
distribution will yield a physically meaningful distribution function. Second, in this analysis 
we have assumed specific forms for the light distribution and the mass profiles obtained will 
be sensitive to these assumed forms. However, as Figure [3] illustrates, our profiles constitute 
reasonable representations of the observed data - but it is worth recalling that the very 
innermost light profiles of the dSphs are often poorly defined. 

Our representation of the velocity dispersion profiles as smooth functions which are flat 
in the innermost regions may mask features which are visible in the profiles at marginal 
significance. It has been known for many years that the stellar populations in dSph are 
complex, with the implications of this complexity evident in many analyses, yet difficult to 
model fully without very large data sets. With respect to a superposition of two populations, 
one should be careful when comparing this to the observed star formation histories. Episodic 
star formation with clearly distinguished periods of star formation has been detected only in 
Carina; all other dSphs with extended star formation histories show evidence of long-lasting 
activity without obvious pauses. So one is not really dealing with, e.g., only two distinct pop- 
ulations (as is sometimes claimed in the literature). It can be shown, however, that younger 
and/or more metal-rich populations are more centrally concentrated (e.g., Harbeck et al. 
(2001)), and that these po pulations also tend to exhibit lower vel o city d ispersions (Tolstoy 



et al. (2004)). Recently, iMcConnachie. Penarrubia. &: Navarro I (120071 ) have recalled the 



general point that a rising dispersion profile (e.g. Leo I) might arise from the superposition 
of two populations in the dSph with different velocity dispersions and spatial scale-lengths. 
While this is certainly an interesting suggestion, we note that all our assumed functional fits 
to the observed dispersion profiles are statistically consistent with the observed data. 

Given these caveats, we conclude that the Jeans analysis demonstrates that the observed 
velocity dispersion profiles and cored light distributions of dSphs are likely to be consistent 
with their inhabiting dark matter haloes with central cores. We note that there is no reason 
that the photometric scale length is exactly the underlying mass scale length. From the 
mass modelling, the lower limit on mass density core size is constrained to be at most a 
factor two smaller than the observed luminosity core in UMi and Fornax. Similarity, rather 
than exact equality, of the two scales is what is relevant here. As discussed below, however, 
cusped mass distributions can also reproduce the observed data on the light profile and 
velocity dispersion profile. To determine the actual slopes of the inner dark matter density 
profiles, further information, either in the form of larger velocity data sets which permit 
full distribution function modelling (see below) or complementary dynamical evidence is 
required. Fortunately, in two special cases, those of Ursa Minor (UMi) and Fornax, there 
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is additional information that enables us to distinquish between shallow and steep internal 
mass density profiles. 

In UMi, an otherwise very simple system fr om an astrophy s ical p erspective, an extremely 
low velocity dispersion sub-structure exists. iKleyna et al. I (120031 ) explain this as a star 
cluster, which has become gravitationally unbound (the normal eventual fate of every star 
cluster), and which now survives as a memory in phase space. Why does it survive in 
configuration space? The group of stars has the same mean velocity as the systemic velocity 
of UMi, so it must orbit clo s e to t he plane of the sky, and hence through the central regions 
of UMi. As IKleyna et al. I (120031 ) show, persistence of the cold structure is possible only 
if the tidal forces from the UMi central mass gradient are weak. In fact, survival of this 
phase-space structure in configuration space requires that UMi has a slowly varying inner 
mass profile, that is, a core, rather than a cusp. 

The Fornax dSph galaxy has five surviving globular clusters. The orbit of a compact 
massive system, such as a globular cluster, should decay due to dynamical friction as it 
orbits through the background dark matter halo particles. The rate of this orbital decay 
is faster in a steep (cusped) dark-matter density profile, and is slower in a shallow (core) 
dark-matter density profil e. The (projected) d istribution of the surviving clusters has been 



(120061 ). who show that a cored mass distribution is 



analysed most recently by lGoerdt et al. 
strongly preferred. 

In summary, while Jeans' equation dynamical analyses cannot be assumption indepen- 
dent, in both cases where some independent information is available, shallow (cored) mass 
distributions are preferred. There is no case where a steep (cusp) distribution is required by 
the data. Applying Occam's razor, we therefore assume for the remainder of this analysis 
that all dSph have similar underlying dark matter mass profiles, which are cored. We adopt 
this specific case- dependent result as a general result since it provides a natural context for a 
characteristic length scale, which is suggested by Figure 1. A single mass model, supported in 
two specific cases with suitable data, links two otherwise disparate results, one photometric, 
one kinematic. 



3.3. King- model dynamical analyses 

Dynamical analysis requires some simple and mathematically smooth functional de- 
scription of the spatial distribution of the stellar tracers of the gravitational field. While 
any smooth function is adequate, those most frequently used include Plummer models and 
King models. Each is convenient, but can mislead if the parameters in the adopted fitting 
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function are (over) interpreted as having physical meaning. 

The King model (King 1966; see also Binney & Tremaine 2000) is physically valid for a 
self-gravitating system with a velocity distribution function which is a lowered Maxwellian, 
i.e. approximates an isothermal distribution at small radii, with an imposed small core - to 
avoid an unphysical divergence - and an imposed cutoff at large radii - to prevent infinite 
extension, and infinite velocities. This model is a good description of a stellar globular cluster, 
in which mass follows light. The equilibrium velocity distribution function which underlies 
this model is that generated by cumulative long-range gravitational interactions between the 
component stars, which brings a system of N stars with rms velocity dispersion v, scale size 
R, and corresponding crossing time t cross into dynamical equilibrium in a characteristic time 
« N/8ln(N) x t cross ;t cross w R/v. It is therefore naturally appropriate to any system in 
which the astrophysical lifetime is long compared to the dynamical relaxation time, and in 
which stars are orbiting in a gravitational potential generated self-consistently from their 
mass. 

We review the relevant range of applicabilities in Table 2. Only the systems in that 
Table with dynamical age > 1 are amenable to a physically meaningful King-model analysis, 
where we define 'dynamical age' as the ratio of the astrophysical age to the relaxation time. 
dSph galaxies have a dynamical age 3-4 orders of magnitude outside this range of validity, 
being too young in a dynamical sense. 

We may ask what would be the stellar mass of a system with the internal velocity 
dispersion and spatial scale size of a representative dSph galaxy, and which has its dynamical 
relaxation time less than the Hubble age of the Universe, so that the physical conditions 
appropriate to establish a King model are in place. The result is that one requires a stellar 
system of 5 x 10 10 stars. For a plausible stellar mass to light ratio, this implies a galaxy of 
luminosity > 1O 1O L . Observed dSph galaxies are many orders of magnitude less luminous. 
Thus, one does not expect that a King model will be a physically valid description of a 
dSph galaxy, although, as Table 2 shows, one does expect such models to be reasonable 
approximations for massive star clusters. While it may still be convenient to use a King 
model as a fitting function for a galaxy, it is invalid to interpret the two free parameters 
of that fit - the 'core radius' and the 'tidal radius' - as meaningful physical properties of 



the galaxy. IWul (120071 ) also shows that King models are inadequate descriptions of available 



dSph kinematics. 

More generally, there are unavoidable consistency requirements in any mass follows 
light model. In any model where mass follows light the projected velocity dispersion must 
be maximum at the centre, and then fall monotonically. For a well-mixed (star cluster) 
system, the velocity dispersion will decrease by roughly a factor of two over three core radii. 
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This is an unavoidable requirement for any mass follows light system, and is observed in star 
clusters (Figure 2). Such a velocity dispersion profile is not required by data in any well- 
studied galaxy, however small, further emphasising the intrinsic difference between (virial, 
King) star clusters and galaxies. 



3.4. Time-dependent kinematics and radial range of valid analysis 



Determination of the mass distribution in the outer parts of a dSph satellite galaxy re- 
mains a complex and data-starved challenge. The (very few) tracer stars at large radii occupy 
the extreme limits of the kinematic distribution function, where simplifying assumptions are 
least reliable. At some radius tidal forces from the Milky Way must become important, 
violating the equilibrium dynamics assumption, though that radius depends on the a priori 
unknown dSph dark matter distribution and the (usually poorly constrained) orbit of the 
dSph around the Galaxy. 

Very many studies are available predicting the effects of time-dependant tides on the 
structure and kinematics observable in the outer parts of galaxies. However, most such stud- 
ies are idealised, and consider o nly single-component models - i.e. either no dark matter, or 

only dark matter, as one prefers ( Johnston. Sigurdsson Sz Hernquist Ill999 ; Johnston. Choi &: Guharthakurl 
20021 ; ISohn et al. II2006I ). This makes any comparison to observations of the stellar density 



distribution or of the stellar kinematics in a dSph galaxy with a dark matter halo at best 
problematic. One recent example, among several, which illustrates th e richness of t he po - 
tential tidal effects on stars orbiting in dark-matter potentials is that of iRead et al. I (120061 ). 
which is one of the few studies to consider external tidal effects on two-compo nent dSph mod- 



els in which stars orbit inside a dark-matter halo. In another study of tides, iKlessen et al. 



( 120031 ) model the available data for the Draco dSph, under the assumption that there is no 
dark matter (ie a single-component model), so that the dSph is unbound. They deduce that 
the present smoothness and small line-of-sight depth of this galaxy make unbound models 
impossible, and conclude that dynamically dominant dark matter is required. 

The observed kinematic and spatial distributions in the outer parts of the dSph remain 
poorly determined by observations. As may be seen in Figure 2 , there is some evidence for 
low-dispersion (cold) outer populations ( Wilkinson et al. I ( 2004)), and als o for fl at velocity 



dispersion profiles to large distances ( Sohn et al. (12006 ) ; Munoz et al. ( 2005 )). both of 



which are inconsistent with simpl e tidal disrup t ion ef fects, particularly since most dSphs show 



no evidence of apparent rotation iKoch et al. I (120061 ) , another p rediction of ti d al disru ption 
models (Read et al. 2006). Carina may be an exception, as iMunoz et al.l (j2006al ) have 



recently detected a velocity gradient on large scales (beyond the nominal King-model fit 
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'tidal' radius) in this dSph. Photometric studies continue to be inconclusive, with some 



suggesting a char acteristic signature of tidal distortions (e.g. Ilrwin fc Hatzidimitriou 111995 



Sohn et al. 1120061 ). but with later studies of the same galaxy failing to find any signal from 
improved data (Odenkirchen et al (2001), Segall et al (2006)). Considerable uncertainty 
about the dynamical state of the three dSphs nearest to the Galactic centre however remains, 
with the closest one - Sgr - being manifestly disrupted. 

Fortunately for our present purposes, which are concerned with the properties of the 
mass distribution in the inner regions of the dSph, the dynamical state of the very outer 
parts of the dSph is relatively unimportant. The complexities noted do mean any discussion 
of total masses is currently impracticable. 

We also note for now that if indeed it is shown from future studies that dSph galaxies 
are not strongly dark matter dominated, but are star-cluster like stellar systems whose 
structural and kinematic properties have been inflated in some way, the conclusions of this 
paper concerning a minimum scale length on which dark matter is seen to concentrate are 
inevitably and drastically strengthened. The minimum dark matter clustering scale would 
have to extend beyond that derived here, and become of order lkpc. 



3.5. Distribution function modelling 

A known limitation of Jeans' equation moment analyses is that (at least) some of the 
dSph galaxies show complex stellar populations, so that adoption of a single dispersion profile 
and single length scale is necessarily a simplification. In all cases these Jeans' moment 
analyses are applicable only over the range where simple functions provide an adequate 
description of the underlying galaxy, which in general is limited to one to two physical 
(luminous) scale lengths. More complex behaviour, preventing valid application of such 
simple models, is seen in the very outer parts of most dSph studied to date, where also 
identification of member stars becomes increasingly uncertain. In general, a more robust 
analysis requires significantly more information. 

As larger data sets become available distribution function modelling can supersede use of 
the Jeans' moment equations. In distribution function analyses one proceeds by constructing 
parameterised equilibrium dynamical models, allowing the dark halo shape and mass, and the 
tracer velocity anisotropy to vary. From these models one can determine model distributions 
of the observable line of sight velocity. These models are then convolved with observational 
errors and an orbital velocity distribution for binary star systems appropriate to a dataset 
of interest, to predict observable velocity distributions at every point across the projected 
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galaxy. It is then straightforward to determine the best fitting models using the individual 
stellar velocities, without the need to degrade the data into moments (dispersions). 



Models of this type have been applied by iKleyna et al. I (120011 ) and by I Walker et al 
( 120061 ). and are being developed further for future application to all large available data 



sets. 'Mass follows li ght' models for Draco were ruled out at the 2.5er confidence level using 



this type of analysis (jKleyna et al. I (]200ll )). Constant anisotropy mode l s of this type favou r 
rather shallow halo inner mass profiles with p oc r -0 5 



Maeorrian 


2003; 


Koch et al. 


2006) 



A form of distribution function modelling has recently been used by lPenarrubia. McConnachie fc Navai 



(120071 ). who apply a methodology developed by Lokas (Lokas 2002; Lokas, Mamon & Prada 
2005) to fit observed dispersion profile data adopting a King-model for the luminous galaxy, 
and embed this inside an assumed NFW dark halo. This analysis requires, as is usual in such 
fitting of NFW models, very considerable dark masses associated with the observed dSph. 
NFW virial masses for the dwarfs considered in this paper are all larger than 10 9 M , with 
that for Draco having log(mass)=9.8 M . This form of modelling breaks the degeneracy in 
the mass determination discussed above essentially by requiring that the assumed dark halo 
be of NFW form. 



Wul (120071 ) has recently applied more general distribution function models to a rebinned 
version of the data for Draco and UMi shown in our Figure 2 here (excluding our outer- 
most data), and for other published data for the Fornax dSph. Wu's analysis considers 
both two- and three-integral distribution-function models, assuming the galaxies are stable 
axisymmetric systems embedded in spherical dark-matter potentials. He includes a range of 
forms for the underlying gravitational potential, considering constant density models, con- 
stant mass-to-light ratio models, dominant central point-mass (massive black hole) models, 
NFW models, power-law and isochrone models. Wu's isochrone model is closest in form to 
those derived here in our Jeans' analysis (Figure 4), having an inner cored mass profile out 
to some radius determined by fitting to the data, beyond which the profile of the dark matter 
distribution steepens. 

Wu's analysis shows that all of the NFW, power-law and isochrone models are consistent 
with the data, while the constant density, constant mass-to-light ratio and dominant central 
black hole models are strongly ruled out. This conclusion is in very good agreement with 
results of our analysis in this paper, where we use additional information to prefer cored 
mass models (isochrone- like potentials) over the NFW and pure power-law cases. Wu's 
characteristic scale lengths at which the underlying dark matter density breaks below the 
cored inner distribution are 500pc, 200pc and 900pc for Draco, UMi and Fornax, respectively. 
These results are quite consistent with the results of the simpler Jeans' analysis we present 
here. 
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3.6. Mass distributions in dSph galaxies 

There are several Jeans' equation analyses of dSph kinematic data which are fully de- 
scribed in the recent literature. Fig ure 4 summarises the r esults of these Jeans equation 



models for several of the dSph from I Wilkinson et al. I (120061 ). with more recent results for 



Leo I from Koch et al (2006) and for Leo II from Koch et al (2007), with in each case 
the simplest possible assumptions for the velocity distribution, namely that it is isotropic 
at all radii. It is apparent that the models are invalid at large radii, where an unphysical 
oscillation in some of the mass profiles is derived. In the inner regions the models are well- 
behaved, and reproduce the overall shape of the observed dispersion and light profiles. As 
illustrated in detail in, for example, Koch et al (2006; especially their figure 11), for Leo I, 
both cored and cusped mass models can provide acceptable agreement with the data for a 
suitable value of a constant anisotropy, and excellent agreement when allowing a radially 
variable stellar orbital anisotropy. As described above, although the Jeans models alone 
cannot distinguish between cored and cusped mass distributions, in those two cases where 
additional information is available cored density profiles are preferred. 

In Table 3 we summarise the mass determination results available. The most robust 
numbers are the cumulative mass within the extent of the kinematic data, and the associated 
mean density and outer circular speed (columns 5, 6, 8 in Table 3). The central density is 
more model-dependent. As an illustration of the range of likely values, column 7 gives the 
mean density within lOpc in the case in which the dark matter has a density profile which 
goes as p oc r~ l throughout the volume occupied by the stellar distribution. If the dark 
matter has approximately constant density out to some break radius, as our models prefer, 
then its central density will of course, be comparable to the mean density quoted in column 
(6). Table 3 also provides an estimate of the mean phase space density within the half-light 
radius, which is defined as 3/(8irGr^a). This estimate assumes that the half-mass radius is 
comparable to the observed half-light radius and that the stellar velocity dispersion is related 
to the total mass through a 2 = GM/r^. Since neither assumption is strictly valid (columns 
(3) and (5) show that this yields an underestimate of the true mass) column (9) should be 
interpreted as an order of magnitude estimate only. 

The total masses within the optical radii of the dSph galaxies have been suspected for 
some years of showing a remarkably small range. Mateo (1998) showed that the available 
data at the time was consistent with an apparent minimum dark halo mass, within the 
optical galaxy, of order 1O 7 M . This relationship was extended and developed by Cote et al 
(1999), who showed it also applied to availab le data from M3 1 satel lites. Later updates were 



provided by I Wilkinson et al. 



(120061 ) and by iGilmore et al. I (120061 ) . The current version of 



that relationship is shown in Figure 5. Remarkably, the Mateo proposal has survived an 
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increase in the dynamic range of the sample by an order of magnitude in both axes, and has 
become better-established, and of lower scatter, as newer data have become available. 

The Mateo plot presents total dark masses within the optical radii. A total mass is the 
integral over a scale length, a mass density profile, and a central mass-density normalisation. 
Each of those parameters is addressed in this paper. The photometric scale lengths are 
summarised in Figure 1, which illustrates that dSph galaxies have both a minimum scale 
size (~ 120pc), and a rather small range of scale sizes. The available kinematic data and 
their analyses are discussed above, summarised in Figure 4. That shows that each galaxy 
studied, under the assumption of a valid Jeans'-equation analysis, has a similar mass profile, 
both in shape and in normalisation: in at least two cases it is probably cored, and the 
central density, assuming a cored profile, is then very similar for all cases analysed to date. 
For clarity, we note that have no robust independent proof that the photometric scale length 
is exactly the underlying mass scale length. From the mass modelling, the lower limit on 
mass density core size is constrained to be at most a factor two smaller than the observed 
luminosity core in UMi and Fornax. Similarity, rather than exact equality, of the two scales 
is what we are assuming here. Constancy of the scale length, the normalisation, and the 
density profile parameters naturally explains the Mateo plot, which is their product. 



3.7. Are these derived properties robust? 

The photometric (length-scale) data are described above, and shown to be robust for 
star clusters, and conservative relative to the present conclusions for galaxies. 

There are two dynamical effects of relevance to bear in mind. In dense star clusters 
the very steep and highly time-dependent gravitational potential gradient between binary 
stars (especially close binaries) provides a force which could eject dark matter particles. 
Essentially, close binary stars will orbit through the cluster threshing the local potential 
gradient, clearing a dynamical tunnel through phase-space, and ejecting any dark matter 
particles which might have been present. Such an effect will be particularly important in the 
central regions of star clusters, which are continually occupied by close binaries. This effect 
is worth additional consideration, and suggests that more extensive numerical modelling 
and kinematic studies of the outer parts of diffuse star clusters is worth while. Significant 
dynamical evolution will be irrelevant in systems whose internal dynamical relaxation time 
is much longer than its age. We show above that all the small galaxies of interest here in 
fact will be immune to this effect. 

The second important consideration concerns the validity of the assumed steady-state 
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dynamical analysis. There is continuing debate and study of the possible relevance of time- 
dependence in dynamical analyses of dSph galaxy stellar kinematics. Time-dependent tidal 
disruption is manifestly dominant in the most nearby dSph galaxy, Sgr, and is probably 
relevant for other dSph galaxies within a few tens of kpc from the Galactic centre (ComaBer 
and UMall are the obvious candidates). In the more distant galaxies there is no kinematic 
evidence that the internal stellar kinematics in their central regions are in any way affected 
by external Galactic tides (IKoch et al. M2006I ) . Most - but not all - kinematic studies remain 
limited to moderate sample sizes, so that necessarily statistical and model-dependent re 



moval of interlopers and outliers can affect conclusions somewhat (eg I Klimentowski et al 



( 120071 ) ) . Similar conclusions have very recently been derived by Wu (IWu l 2007h m his 
parametric re-analysis of published data for Draco, Fornax and UMi. Wu concludes: " 
Because [his distribution function] models can fit both radial velocity profiles and surface 
number density profil es, the so-called "extra-t i dal ex tensi ons" in the surfac e numb er den- 
sity profiles found by llrwin fc Hatzidimitriou I (119951 ) and I Wilkinson et al. I (120041 ) do not 



require any special ad hoc explanation. Thus it is not valic 
dence of tidal stripping, as proposed by some authors (e.g., 


to consider them as the evi- 


Martmez-Delgado et al. 


2001; 


Gomez-Flechoso & Martmez-Deleado 


2003; 


Munoz et al. 


20051)." 



We note for completeness that there are continuing efforts to apply pure tidal models 
- i.e. there is no dark matter, and the velocity dispersions are inflated by tides - to dSph 
kinematics. If these models can be proven to be relevant, the implications for the nature of 
dark matter are quite profound. Substantially longer minimum length scales and substan- 
tially lower maximum mass densities even than those derived here will be required, providing 
quite extreme constraints on the nature of dark matter. It remains far from clear that such 
mass-follows- light models are consistent with the robust evidence on cosmological scales that 
dark matter dominates the Universe. 



4. Discussion and Implications 



There is strong, and strengthening, observational evidence for four conclusions concern- 
ing the smallest galaxies and typical or brighter (globular, nuclear) star clusters. These 
two families of objects co-exist over the range of stellar absolute magnitudes from — 15 ^ 
My < —4, corresponding to luminosities from ~ 10 9 L m > L > 10 3 L M , with the upper limit 



corres p onding to the brightest known star clusters ISeth. Dalcanton. Hodge fc Debattista 



(bood ) ; Irlilker etal I (bo07h : lEvstigneeva. Gregg. Drinkwater. fc Hilkerl (l2007h ]. the lower limit 
to the least luminous galaxies as yet discovered (IBelokurov et ajj (120071 )). The number of 
well-studied objects known within these limits in both classes has increased substantially 



in recent years, pr oviding an adequate sample to identify s ystematics, and t o test earlier, 
provisional trends (jCote. Mateo. Olszewski fe Cook I (119991 ); I Wilkinson et al. I (120061 )). 



We conclude the following: 

ONE: Over this substantial dynamic range, there is a clear bimodality in the size 
distributions of the two families of object: all smaller objects are star clusters, and have 
characteristic scale sizes ^ 30pc (< 10 18 m), while all larger objects are galaxies, and have 
scale sizes for their luminous (stellar) components > 120pc (> 4.10 18 m). 

TWO: Where kinematic studies exist [essentially in the central luminous regions], there 
is a clear distinction between the phase-space distribution functions of star clusters and galax- 
ies. At a given (stellar, baryonic) luminosity, galaxies have phase-space densities typically 
two to three orders of magnitude lower than do star clusters. In all adequately-studied cases, 
the galaxy's luminous stellar component is embedded in a more extended dark matter halo. 
Star clusters over the whole mass range are well- described by the virial theorem, and show 
no evidence for any (dynamically significant, extended) dark matter halos. 

THREE: In the two specific galaxies where both detailed dynamical analyses are fea- 
sible, ie, substantial kinematic data across the face of the galaxies is available, and where 
independent evidence to break the core/cusp/velocity- anisotropy mass degeneracy exists, 
the derived dark mass distribution has a shallow density profile. For a density distribution 
p describable as a function of radius r as p(r) oc r~ a , the data imply that the power-law 
index a < 0.5, and with a being consistent with zero in the innermost regions. Simplicity 
argues for this being the general case. Under that assumption of generality, the minimum 
photometric length scale can be interpreted as comparable to a minimum length scale for the 
clustering of dark matter. The lower limit on the size of a mass density core is constrained 
to be at most a factor two smaller than the observed luminosity scale, while our dynamical 
modelling suggests it is not much, if at all, larger. 

FOUR: Using the mass-model assumed in point THREE, the derived mean mass 
density of dark matter within one to two half-light radii for all the galaxies is pdm ~ 
5GeV/c 2 cm -3 (about 0.1 M pc~ 3 ). If we do not adopt the results of point three, but 
allow a cusped mass model, the derived maximum mass density of dark matter within lOpc 
of the centre is p ma x,DM ~ 2TeV/c 2 cm -3 (about 60 M pc~ 3 ). 

The combination of the small range of observed scale sizes, together with an apparently 
standard (isochrone-like) form for the derived dark mass density profile and its normalisation, 
naturally explains the observed relation that all dSph galaxies have similar total dark mass 
within their optical radii. 
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These systematic properties have obvious implications for the nature of dark matter and 
galaxy formation in small halos, which we note very briefly here. The simplest case is that our 
cored mass profiles are the unmodified outcome of formation of small halos. An alternative 
is that dark matter profiles on these small scales have been drastically restructured by some 
astrophysical (feedback) process. We very briefly note each in turn. 

4.1. Implications for ACDM galaxy formation models 

One immediately asks if our observed length scale and characteristic upper mass density 
derived here for the intrinsic nature of dark matter is consistent with the extremely high- 
quality agreement between ACDM models and large-scale structure, assuming the shallow 
small-scale profiles are an intrinsic feature of dark matter halo formation. Few (if any) cos- 
mological observations of large-scale structure have resolution on small enough scales to be 
sensitive to the truncation of the small-scale power spectrum on scales of order lOOpc implied 
here, so there is no disagreement. Galaxy formation models inside the ACDM paradigm how- 
ever have considerable difficulties matching observations on small scales. The well-known 
'satellite problem' is an example, as is the 'cores vs cusps' debate. The amount of structure 
on sub-kpc scales predicted by simulations will be very drastically modified by our conclu- 
sions here, which imply that the primordial power spectrum is truncated at small physical 
length scales. In both the satellite counts and the core-cusp debates, the discordance be- 
tween data and numerical models will be greatly reduced in simulations which include the 
minimum scale cutoff suggested by our results. Numerical studies are underway to quantify 
this. 

A fundamental astrophysical, rather than dark matter, puzzle which is not simply ex- 
plained is the now well-established bimodal size distribution function illustrated in our Fig- 
ure 1. Star clusters have a maximum half-light scale size of some 30pc, while the charac- 
teristic minimum scale length associated with the stellar systems which occupy dark matter 
halos is some 4-5 times larger. The dark matter mass profiles derived here have a scale 
length always a factor of several longer than the scale on which self-gravitating star clusters 
form. But why do stellar systems not form with all possible scale lengths inside the shallow 
dark matter potentials which we observe kinematically? It is a generic prediction of ACDM 
galaxy formation models that the baryonic component should cool and collapse by a factor 
A -1 PS 10 more than does the (presumed self-interaction-free) dark matter. Why is this not 
seen in these very low-mass halos? 

It seems that stellar populations in these shallow potentials expand to a scale size 
comparable to that of the underlying potential. An attempted explanation is beyond our 
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present purposes. We note however th a t there are many studie s of feedback in shallow poten- 



tial w el ls (see eg. iDekel fc Silk! (119 86); Sil k. Wvse. fc Shields (119871 ): iRead. Pontzen fc Vie! 



( 120061 ); iMayer. Katzantzidis. Mastropietro. &; Wadsleyl (120071 )) which discuss the clear im- 
portance of energy balance in driving gas loss an d in expanding extant s tellar systems. We 
have attempted an explanation of this ourselves (IRead fc Gilmord (120051 )). concluding that 
plausible astrophysical feedback capable of leaving a dSph-like remnant cannot convert a 
cusped DM profile into a cored DM profile, but extended shallow (exponenti al) dSph lu- 
minosity can be readily formed. More generally IRead. Pontzen fc Viell (120061 ) provide an 
extensive discussion and introduction to the literature. Several of these models remove all 
residual gas in an early star formation episode, with the lost gas-mass reducing the binding 
energy of the residual stellar system, allowing it to expand into the dark matter potential. 
A further general constraint on these promising models is that many dSph galaxies have 
supported extended star formation, and some show evidence for internal chemical element 
self-enrichment, sometimes over most of the age of the Universe. Derived star formation 
rates tend to be very low - of order one star per 10 5 years - so that single disruptive gas 



evacuation events are not part of solution space (e.g.. iHernandez. Gilmore fc Valls-Gabaud 



2000l ; ICarigi. Hernandez fc Gilmore 1 120021 ) . Such disruptive events are of course consistent 



with the apparent absence of stellar systems occupying dark halos of lower mass than those 
we observe, should they be more prevalent on smaller mass scales that we study here. 



4.2. Implications for the nature of dark matter 



An adequate discussion is beyond our intention in this paper, which is to establish the 
observational evidence. The local dark matter mass density near the Sun was determined us- 
ing di stribution function modelling to have the value pdm ~ 0.3GeV/c 2 cm~ 3 (about 0.01 M 



pc 3 ) (jKuijken &; Gilmore Ill989l . Il99ll ). Our maximum mass density and minimum physical 



length scale derived here imply particle number densities, and self-interaction cross sections 
(if one assumes self-interaction is the physical cause of the present-day length scale) which 
are readily calculable for any specific particle class. We note only the obvious conclusion 
that the very low maximum mass density derived here is challenging for models of dark 
matter which are dominated by massive (of order TeV) particles, such as those predicted by 
some supersymmetric theories. Such particles would be required to have a spatial number 
density < 1 cm -3 even if the mass density profile is cusped, and orders of magnitude lower 
if our cored profiles are appropriate. Much lower-mass (/ieV) particles, such as the axion, 
would have correspondingly higher number densities. Rathe r interestingly, intermediate - 
mass (keV) sterile neutrino particle s have been discus s ed (see eglDodelson fc Widrow I (119941) ; 
Abazaiian. Fuller, fc PateTT (hoOlh : iKusenko I kood \: iBiermann fc Munvanezal (bo07al lbhks 
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relevant in just the spatial and density range we have derived here. 
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Table 1. Observed Properties of the Established Milky Way dSph Satellites and 
Candidates. Quoted half-light radii are derived from Plummer model fits from the 
identified source, except for those cases identified as core radii. Core radii, from King 
model fits, are lower limits on the half-light radius. 



C eject 


r 


Galacto-centric 


Half-light Radius 


name 




distance (kpc) 


pc 


Sgr a 


> 2 x 10 7 


24 


> 500 


Fornax 1 


1.5 x 10 7 


140 


400 (core) 


Leo I 2 


4.8 x 10 6 


240 


330 


Sculptor 3 


2.2 x 10 6 


80 


160 (core) 


Leo II 4 


7 x 10 5 


230 


185 


Sextans 5 


5.0 x 10 5 


85 


630 


Carina 6 


4.3 x 10 5 


100 


290 


UrsaMinor 7 


3.0 x 10 5 


65 


300 (core) 


Draco 8 


2.6 x 10 5 


80 


230 


CVn I 9 


> 1 x 10 5 


220 


550 


Hercules 10 


> 2 x 10 4 


140 


310 


Bootes 11 


> 2 x 10 4 


60 


230 


Leo IV 10 


> 1 x 10 4 


160 


150 


UMa I 12 


> 1 x 10 4 


100 


290 


CVn II 10 


> 7 x 10 3 


150 


135 


UMa II 13a 


> 3 x 10 3 


30 


~125 


ComaBer 10 


> 2 x 10 3 


45 


70 


Segue I 1[)b 


~ 10 3 


25 


30 


Willman I 14b 


~ 10 3 


40 


20 



a Sgr certainly and UMa II probably, the objects closest to 
the Galactic centre, a r e associated with ext ended tidal streams: 
iBelokurov et all l l2006ah : Fellhauer et"aH [20071. The quoted param- 
eters in both cases are highly uncertain. 

b Nature uncertain: suspected to be a globular star cluster 

^Walker et al. [ll2006ll ; (21 iKoch et al. I ll2006h: (31 iMateo I lll998ll; 

Koch I 
200C); (7) 
. -20041): (91 

I ll2006a): llbata et al. I I2OO6II; (101 Belokurov et all 
Belokurov et al. I 2006 bll :^unoz et al. T ~i 2006bll: (12) 
lyyillman et al. I ll2005all: iKlevna et al. I ll2005ll ; (131 Zucker et alH 
ll2006dl : lGrillmairN2006ll : (14l lWillman et ahl ll2005bli 



IWestfall et al. I COOrj), 
d2007f); (51 IKlevna et^O j20o3 T 




Irwin fc Hat zidimitriou N1995II: (8) IWilkinson et al 



Zuck er et al. 
J2007h; (11) 
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Table 2. Dynamical relaxation times for dSph and star cluster densities 



Star 


Radial 


1-D velocity 


Crossing 


Relaxation 


dynamical 


object 


numbers 


scale (pc) 


dispersion (km/s) 


time (yr) 


time (yr) 


age 


type 


100 


2 


0.5 


4 X 10 6 


10 7 


£ i 


open cluster 


10 5 


4 


10 


4 X 10 5 


4.10 8 




median glob cluster 


10 8 


10 


30 


3 x 10 4 


2.10 10 


~ 1 


nuclear/large glob cluster 


10 12 


10 4 


300 


3 x 10 7 


10 17 


io- 7 


gE gal 


10 9 


10 3 


50 


2 x 10 7 


10 14 


io~ 4 


dE gal 


10 8 


400 


10 


4 x 10 7 


10 13 


10- 3 


dSph 
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Table 3. Observed velocity dispersions for Milky Way satellites, and derived masses and 
densities. The columns are: (1) Galaxy name; (2) Total velocity dispersion from extant 
data (as quoted in associated reference); (3) Crude mass estimate based on total velocity 
dispersion; (4) Radial extent of mass models - the radial extent quoted is either the actual 
region within which the mass has been calculated (as given in the associated reference) or, 
in cases where this radius is unavailable, the nominal King limiting radius (denoted rn m ); 
(5) Mass within radius (4) based on modelling of extended dispersion profile; (6) Mean 
density within the radius in column (4); (7) Mean density within lOpc assuming mass is 
distributed as a power-law with p oc r _1 ; (8) Circular speed at edge of data = \/GM^^Jr\ 
(9) Estimated mean phase space density within half-light radius = 3/(8TcGr^a). See text 

for an explanation. 



galaxy 






fraax 


M DM 


PDMffmax) 


PDM,cusp(10pc) 


^circ C^*max ) 


Phase-Space Density 


name 


km s -1 


Mq 


kpc 


Mq 


GeV/c 2 cm" 3 


GeV/c 2 cm" 3 


kms -1 


Mokpc _3 kms _3 


Fornax 1 


11.1 ± 0.6 


1 x 10 7 


1.5 


~ 3 x 10 s 


~ 0.7 


1.1 x 10 3 


29 


2 X 10 4 


Leo I 2 


9.9 ±1.5 


8 x 10 6 


r lim ~ 0.9 


3 - 8 x 10 7 


~ 0.3 - 0.9 


310 - 830 


12-20 


3 x 10 4 


Sculptor 3 


7-11 


2 x 10 6 


rilm ~ 1-8 


£ io 7 


~ 0.015 


27 




2 x 10 s 


Leo II 4 


6.8 ±0.7 


2 x 10 6 


riim ~ 0-5 


3 x 10 7 


~ 1.8 


~ 930 


16 


1 x 10 s 


Sextans 5 


8 


9 x 10 6 


0.8 


~ 3 x 10 7 


~ 0.5 


390 


13 


9 x 10 3 


Carina 6 


7.5 


4 x 10 6 


Him ~ 0.8 


~ 4 x 10 7 


~ 0.7 


520 


15 


4 x 10 4 


UrsaMinor 7 


12 


1 x 10 7 


0.5 


> 6 x 10 7 


~ 4.5 


2.1 x 10 3 


23 


3 x 10 4 


Draco 7 


13 


9 x 10 6 


0.5 


> 6 x 10 7 


~ 3 


1.7 x 10 3 


22 


4 x 10 4 


Bootes 8 


6.6 ±2.3 


2 x 10 6 




~10 7 








8 x 10 4 


UMa I 9 


9.3 


6 x 10 6 




~ 10 7 








4 x 10 4 



1 Walker et al. I J2006h: (2lKoch et al. 



(6) 
(9) 



Wilkinson et al I d2006h ; (7) Wilkinson et atT ikon 
Willman et al. I ll2005aF > ; iKlevna et al. I (l2o6*5h . 



2(iii(,D; (3) iMateo I jl998h : IWestfall et al. I ll2006h: 

3k 



(4) Coleman. Jordi. Rix. Grebel fc Kochl d2007h; Koch et al. I d2007h : (5) IKlevna et alJfeOQ^; I Wilkinson et al. I d2006h : 



4h: (8) iBelokurov et al 



f na et al,l 
^006rl):f 



Munoz et al. 



I |2006rJ ): 



J i i i i I i i i i I i i i i L_ 

12 3 

log [half-light radius (pc)] 



Fig. 1. — Absolute magnitude M v vs (logarithmic) half-light radius for well-studied stellar 
systems. The filled symbols are objects classed as galaxies, the open symbols and asterisks 
objects classed as star clusters of various types. Red colours indicate objects associated 
with the Milky Way Galaxy, blue colours are objects associated with M31 and green colours 
indicate more distant objects. Red filled triangles are the well-known dSph, red filled circles 
are those recently discovered, with in each case the photometry listed in Table 1 being 
adopted (references are given in the notes). The least-luminous M31 dSph (blue pentagons) 
are from Martin etal (2006), and have ~ 50% uncertainties. Ringed circles highlight the 
probable star clusters Segue 1 and Willman 1, and the object ComaBer. Open red circles 
are Milky Way globular clusters, from the compilation of Harris (1996), except the two 
largest Galactic globular clusters (Pal 5 and Pal 14) which use the most recent data from 
Hilker (2006). The largest globular clusters in M31 are shown as blue open squares, with 
data from Mackey et al. (2006). Green pentagons are globular clusters in NGC 5128 (the 
peculiar elliptical galaxy Cen A; Harris et al. 2002, 2006; Gomez et al. 2006). Green crosses 
represent nuclear star clusters in a range of external galaxies (Bastian et al. 2006), open green 
triangles are young massive star clusters (Bastian et al. 2006; Seth et al. 2006). Asterisks 
are Ultra Compact Dwarfs (UCD) in the Fornax cluster of galaxies (De Propris et al. 2005; 
Mieske et al. 2002; Drinkwater et al. 2003), and in the Virgo cluster (Hasegan et al. 2005). 
For UCD3* in Fornax we adopt the most recent core measurement (22pc) by Drinkwater et 
al. 2003). Not shown individually are the "Faint Fluffy" star clusters found in the disks of 
lenticular galaxies (Brodie & Larsen 2002) which have absolute magnitudes My ~ —7, and 
sizes in the range ten to twenty pc (1.0-1.3 in log(r/J). Sgr is not shown. Half-light size 
definitions and determinations are discussed further in the text. 





Fig. 2. — Observed line-of-sight velocity dispersion profiles for six dSph galaxies. Also shown 
(lower right) is the model predicted dispersion profile for a Plummer model in which mass 
follows light. The lower left pa nel shows the observed velocity dispersion profile for the 
globular cluster Omega Cen from ISeitzerl (119831 ). The similarity between the Plummer 'mass 
follows light' model and the data for Omega Cen is apparent, with a monotonic decrease 
in dispersion from a central maximum. In contrast, the dSph galaxies do not have their 
maximum dispersion value at the centre, and retain relatively high dispersions at large radii, 
indicating extended (dark) mass distributions. 



-33- 




Fig. 3. — Functional fits to the surface brightness profile (top) and velocity dispersion profile 
(bottom) of the Draco (left panels) and Carina (right panels) dSphs used to derive mass 
profiles based on Jeans equations. Similar fits are used for the remaining four dSphs presented 
in Figure 4. 
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Fig. 4. — Derived inner mass distributions from isotropic Jeans' equation analyses for six 
dSph galaxies. The modelling is reliable in each case out to radii of log (r)kpc~ 0.5. The 
unphysical behaviour at larger radii is explained in the text. The general similarity of the 
inner mass profiles is striking, as is their shallow profile, and their similar central mass 
densities. Also shown is an r _1 density profile, predicted by many CDM numerical simula- 
tions (eg Navarro, Frenk & W hite 1997). The individu al dynami cal analyses are describ ed 



in full as follows: Ursa Minor (I Wilkinson et al. 



LeoII flKoch et al. I (120071 )): Leol flKoch et al I f 



and Wilkinson et al in preparation); Sextans (IKleyna et al.l (120041 ) ). 



2004m : Draco flWilkinson et al. I t004\) : 



2006hV: C a rina flwilkinson et al. I hood ) 
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Fig. 5. — An updated Mateo plot. Mass-to-light ratios are plotted versus absolute mag- 
nitudes for Local Group dwarf galaxies, following a style suggested by Mateo [Mateo etal 
(1993);(1998, his fig 9, lower panel)]. The solid line is the relation for a constant mass (dark) 
halo. The modern data shown here extend the original relation by three magnitudes in 
luminosity and an order of magnitude in mass-to-light ratio, while reducing the scatter by 
an order of magnitude. Data are from the tables in the text. Values for Scl, Andll, AndIX, 
UMa and Boo are based on small kinematic samples, and are less certain than are the results 
for the other galaxies. We explain this correlation as a consequence of the characteristic min- 
imum galaxy scale size shown in Figure 1 convolved with the narrow range of mass profiles 
and mean dark matter densities shown in Figure 4. 



